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Abstract: Effective delay of turbulent boundary layer separation could be achieved via 

closed-loop control. Constructing such a system requires that sensor data be 
processed, real-time, and fed into the controller to determine the output. 
Current methods for detection of turbulent boundary layer separation are 
lacking the capability of localized, fast and reliable identification of the 
boundary layer state. A method is proposed for short-time FFT processing of 
time series, measured by hot-film sensors, with the purpose of identifying the 
alternation of the balance between small and large scales as the boundary layer 
separates, favoring the large scales. The method has been validated by 
comparison to other criteria of separation detection and over a range of 
baseline and controlled flow conditions on a “simplified” high-lift system, 
incorporating active flow control. 
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1. INTRODUCTION 

The identification of turbulent boundary layer separation and its control 
are of great practical importance [1-3]. Most applications require some sort 
of closed-loop control to perform efficiently, and should use knowledge of 
the boundary layer “state” (i.e. attached or separated and intermediate states 
of different degrees of complexity) over the entire separation prone region. 
Moreover, the knowledge as to how close the boundary layer is to separating 
should be updated at the highest possible rate to provide feed-back to the 
controller. Several flow features can be used as indicators of the “state” 
including, vanishing skin friction [1], local minimum of the wall pressure 
fluctuations [4-5] and vanishing pressure gradient [6]. The measurements 
techniques often used include cross-spectral information from neighboring 
hot-film sensors [7-12], information from the newly developed wall mounted 
miniature LDV system, and floating MEMS devices [3]. 

An appealing technique is the surface mounted hot-film, due to its 
simplicity, availability and low cost. Transition detection using the standard 
deviation alone or with either the skewness or the flatness of un-calibrated 
hot-film sensors is well established [7-13], as is the identification of 
stagnation points and laminar separation bubble using cross sensor 
information for identifying the “phase reversal” feature [8-13]. 

In a turbulent boundary layer, the phase reversal technique cannot be 
used for reliably determining where separation occurs. Several researchers 
[11,13] have presented data using the correlation coefficient of adjacent hot- 
film sensors for detecting separation. None of the above-mentioned methods 
were successful at determining the turbulent boundary layer separation 
location using the present data set. Furthermore, all methods rely on cross 
sensor information transfer, slowing down the loop. Therefore, the search for 
a new turbulent boundary layer separation detection criterion is warranted. 

The subsequent sections of the paper will describe the experimental 
setup, the proposed criterion, the algorithm developed, and the application of 
the method to several data sets acquired on different configurations and at 
different flow conditions. 


2. EXPERIMENTAL SET-UP 
2.1 The wind tunnel 

Tests were conducted in the Basic Aerodynamic Research Tunnel 
(BART), at NASA Langley Research Center. The BART facility is a low- 
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speed open-circuit wind tunnel, with a 0.71 m high by 1.02 m wide by 3.05 
m long test section. The maximum speed of the tunnel is 60 m/s 
(Re/m=0. 345x1 0 6 ). The wind tunnel is well instrumented and allows optical 
access for Particle Image Velocimetry (PIV) measurements. 

2.2 Simplified High-Lift Model 

The simplified high-lift version of the NASA EET (Energy Efficient 
Transport) airfoil model [14] was designed in a modular manner so that 
zero-net-mass-flux actuators could replace solid regions in the model near 
the leading and trailing edge flap shoulders (Fig. la). The 406.4 mm chord 
model has a 15% chord leading edge flap (LEF) that can be deflected from 
0° to -30° and a 25% chord trailing edge flap (TEF) that can be deflected 
from 0° to 60°. Angle of attack settings for the airfoil, the TEF and the LEF 
were all automated and computer controlled. The model has 78 streamwise 
static pressure taps located at mid-span and two rows of 18 spanwise static 
pressure taps spaced 50.8 mm apart located at x/c=0.35 and x/c=0.94. In 
addition to the static pressure taps, there are nine unsteady pressure 
transducers on the model surface (Fig. lb) and at least one unsteady pressure 
transducer embedded in each actuator cavity for monitoring the pressure 
fluctuations produced by the actuator and correlating the wind tunnel 
experiment with the bench-top actuator calibration tests. The data presented 
herein has not been corrected for wind tunnel wall interference. 

2.3 Hot-film Arrays 

To aid in determining the locations of transition and separation, 48 hot- 
films were installed on the model. On each element (LEF, TEF, and main) 
of the model there are 16 hot-film sensors approximately 50.8 mm to the 
right of the model centerline (and therefore from the steady pressure tap 
locations, Fig. lb). A 16-channel constant temperature anemometer coupled 
with a switch matrix was used to operate and acquire data from the hot- 
films. The hot-films on each group of 16 elements of the model were 
acquired simultaneously using a 16-bit high-speed analog-to-digital 
converter. The data were low pass filtered at 10 kHz and sampled at 25.6 
kHz. The films were operated at an overheat ratio of 1.2. The hot-film 
sensors, 0.4 mm long each, were etched onto a polyimide sheet and then 
bonded to the model. The polyimide sheet covered one half of the span of 
the model. A 0.1 mm step exists at the juncture between the polyimide sheet 
and the model. Body filler was used to fair the step. The frequency response 
of the hot-films was evaluated to be flat up to 5-7kHz. 
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Figure la: EET airfoil with actuators regions represented by lighter gray-level. 



Figure lb: EET pressure tap and hot-film locations, c=406.4mm. 



2.4 Actuators 

An internal piezoelectric fluidic actuator was used in each of the actuator 
regions shown in Figure la. The actuator slots were all inclined downstream 
at an angle of about 30 deg to the surface. The flap actuator is used for the 
data presented in this paper. This actuator has four alternative excitation 
slots. Only one slot was active during each experiment. The slot located at 
x/c=0.757 with the TE flap un-deflected is discussed in this paper. This slot 
is 0.635 mm wide and has 19 segments that are 0.051 m long. The remaining 
slots were sealed using either a water-soluble filler or 0.051 mm thick 
adhesive Kapton tape. 

2.5 PIV Set-Up 

Digital PIV was used to measure the instantaneous flow fields, phase 
synchronized with the actuators’ cycle. The PIV system includes two 1.3 k 
by 1 k cameras, equipped with 105 mm Macro lens, installed side-by-side. 
The field of view from each camera was overlapped to cover a larger region 
above the airfoil. The magnification of the imaging system was about 9:1 
with the measurement plane about 70 mm wide. The interrogation area at 
each grid point had a 24x24 pixel resolution. This corresponds to about 

1.5 mm square at the measurement plane. A maximum overlap of 50% 
between adjacent interrogation regions was used. Smoke (average size of 
1 micron), introduced upstream of the contraction, via a commercial smoke 
generator, was used for seeding. Dual Nd-Yag lasers were used to illuminate 
a light sheet, placed about 50 mm off the model centerline. The laser pulse 
separation was set at 55 microsecond to cover a free stream velocity of about 
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10 m/s. Time averaged data acquired using approximately 100 image pairs 
for each flow condition are presented. 

2.6 Experimental Uncertainty 

The a’s presented are accurate to within ±0.03°. The LEF and TEF 
deflection angles are accurate to within ±0.25°, Cp (excitation momentum 
coefficient, [19]) is accurate to within 20%, the chord Re is accurate to 
within 3%. PIV velocity uncertainty for the current test condition is about 
3%. 


3. ALGORITHM 

If one examines the wall normal distribution of the turbulent eddies in an 
attached turbulent boundary layer [ 15 ], one finds that the smaller eddies are 
active close to the wall while the large eddies are active across the entire 
boundary layer. As the turbulent boundary layer approaches separation, due 
to an adverse pressure gradient, the level of the turbulent activity increases 
and the relative magnitude between the small and large scales, as interpreted 
by the spectra of the near wall velocity [ 5 ], the wall shear-stress or the wall 
pressure fluctuations [ 4 - 5 , 16 - 18 ], is altered such that the lower frequencies- 
larger eddies contain an increasingly larger fraction of the turbulent kinetic 
energy. This process, however, is not steady [ 16 ], but intermittent. Our 
physical interpretation of these findings is that at incipient detachment , the 
small but still energy containing eddies will stop interacting with the wall in 
an intermittent manner. The largest scales would still interact with the wall 
even when the flow is mostly separated, affecting the attached flow periods, 
and therefore, also reinstating the high frequency content of the measured 
wall mounted sensor signal. Hence, a flow separation criterion was 
developed, that is based on the intermittent disappearance of the high 
frequency content from the spectra of a hot-film sensor signal. This criterion 
does not require cross sensor information exchange and could hopefully be 
used also in complex 3D flows. 

The method is constructed of the following steps: 

1 . Short time FFT (STFFT) of the hot-film voltage time history, 

2. High-pass filtering of the FFT coefficients, above a 
predetermined cross-over frequency (to be defined below), 

3. Summation of the high frequency coefficients to generate the 
“high-frequency standard deviation” (HFSDV) of the signal, 

4. Time averaging of all the HFSDV’ s calculated in step 3, 

5. Normalizing the time dependent HFSDV’s by their mean, 
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6. Calculating the standard deviation of the normalized HFSDV. 

The crossover frequency (f c ) was selected based on interrogation of a 
large number of data sets to be f c /Ue»5, where Ue is the velocity at the 
boundary layer edge, calculated from the pressure distribution. The result of 
the above procedure is independent of the flow conditions due to the auto 
normalization. Time windows of 40ms were currently used, at the current 
low speed conditions, but it can be significantly shorter especially at higher 
speeds. A threshold level of 0.225 was found to be adequate for identifying 
baseline as well as controlled separation over a wide range of configurations 
and Reynolds numbers. 

The addition of large or small-scale coherent structures [19-23], for the 
purpose of controlling the turbulent boundary layer separation, complicates 
the separation identification procedure, by adding coherent motion to the 
complex separating boundary layer spectra. However, a separation of scales 
exists between the small-scale motion, which the current criterion searches 
for its disappearance, and the excitation imposed large-scale structures. The 
method also performed well when high frequency excitation was introduced, 
due to the high dissipation rate of these scales. 


4. RESULTS 


4.1 Baseline and controlled flow 

Due to the limited scope of the current paper, only one flow condition 
will be described in detail. Fig. 2 describes the pressure distributions on the 
baseline and controlled flows. The model is at incidence of a=6°, the leading 
and trailing edge flaps are deflected -25° and 20°, respectively and the chord 
Reynolds number, Re c , is 0.24xl0 6 . 

The baseline flow separates at x/c«0.8, i.e. slightly downstream of the 
deflected flap shoulder, based on the vanishing pressure gradient criterion. 
When high frequency excitation at rather high momentum coefficient, 
Cp«1.0%, was applied from the flap actuator slot, located at x/c=0.77 
(8f=20°), the flow partly reattached to the flap. This actuation resulted in 
delayed separation and increased circulation around the entire airfoil. The lift 
coefficient was increased in this case from 1.65 to 1.82 while the form-drag 
coefficient increased from 0.074 to 0.078. 
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4.2 Application of the algorithm 

The mean wall pressures (Cp) could only provide an indirect indication 
of the location of the mean separation location. A more direct indication of 
the mean separation location can be obtained from analyzing the PIV 
measured streamwise velocities presented in Figures 3a and 3b, for the 
baseline and controlled flows, respectively. Note the seam between the two 
partly overlapping images acquired simultaneously by the two PIV cameras. 
One should take into account the wall normal averaging of the PIV 
technique, due to the interrogation cell size. The PIV data indicates, using 
the zeroing of the near wall velocity [1] as guidance, that the baseline flow 
separates at x/c=0.78. According to the same criterion, separation is delayed 
to x/c=0.87, using the control input. 

Applying the separation detection criterion based on the hot-film signals 
and using the suggested threshold level of 0.225, as shown in Fig. 4, 
indicates that the baseline flow separates at x/c=0.78 and that separation is 
delayed to x/c=0.87 in the controlled flow. The agreement between the Cp, 
the PIV data and the new separation detection criterion is satisfactory, 
validating the criterion. The criterion was validated also at different Re c and 
LEF and TEF deflections and excitation locations [ 23 ]. 

Sensitivity tests for the sampling rate, STFFT window width, low pass 
filter setting to simulate a lower band-width sensor, and threshold levels 
were conducted and resulted in low sensitivity. 


5. CONCLUSIONS 

A new turbulent boundary layer separation detection method has been 
proposed, formulated and validated using experimental data obtained over 
the NASA-TAU “simplified high-lift system”. 

The method makes use of the changing balance between small and large- 
scale turbulent structures, favoring the large scale, as the boundary layer 
separates. Raw surface mounted hot-film sensor data was acquired and 
analyzed. The method and the separation detection criterion were validated 
using mean pressures and PIV data over a range of Reynolds numbers, 
geometries, and excitation locations. 

Further study is required in order to enable the use of the new method 
and criterion as part of a closed-loop active separation control system, with 
the aim of efficient-distributed control of turbulent boundary layer 
separation. 
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Fig. 2 Pressure distributions of baseline and controlled flow [19]. Dashed box indicates 
PIV interrogation region of Figures 3. 
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Fig. 4 Separation location based on a threshold level of 0.225 for the flow conditions of 

Figures 1-3. 


